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Abstract. Oscillatory critical amplitudes have been repeatedly observed in hierarchical 
models and, in the cases that have been taken into consideration, these oscillations are so 
small to be hardly detectable. Hierarchical models are tightly related to iteration of maps 
and, in fact, very similar phenomena have been repeatedly reported in many fields of math- 
ematics, like combinatorial evaluations and discrete branching processes. It is precisely in 
the context of branching processes with bounded off-spring that T. Harris, in 1948, first set 
forth the possibility that the logarithm of the moment generating function of the rescaled 
population size, in the super-critical regime, does not grow near infinity as a power, but it 
has an oscillatory prefactor. These oscillations have been observed numerically only much 
later and, while the origin is clearly tied to the discrete character of the iteration, the am- 
plitude size is not so well understood. The purpose of this note is to reconsider the issue 
for hierarchical models and in what is arguably the most elementary setting - the pinning 
model - that actually just boils down to iteration of polynomial maps (and, notably, qua- 
dratic maps). In this note we show that the oscillatory critical amplitude for pinning models 
and the oscillating pre factor connected to the Harris random variable coincide. Moreover 
we make explicit the link between these oscillatory functions and the geometry of the Julia 
set of the map, making thus rigorous and quantitative some ideas set forth in [10] . 
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1.1. Hierarchical models, Branching processes, and the main issue. It has been 
pointed out by several authors for a long time, see e.g. [19], that the discrete character of the 
renormalization group transformation may in principle give rise to a periodic modulation of 
the critical amplitude. This oscillation has been then observed in hierarchical models, see for 
example [10] that deal in particular with Ising and Potts models on diamond lattices. 
The modulation of the critical amplitude in these models turns out to be very small and the 
nature of this phenomenon has not been fully elucidated. Here we consider this phenomenon 
for a particular hierarchical model, the wetting or pinning hierarchical model |9j, because it is 
arguably the easiest set-up and it makes a direct contact with a vast mathematical literature: 
iteration of polynomial maps. 

In fact the partition function of a hierarchical model with pinning potential h G R and 
volume size d n if just / n (exp(/i)), where f n = f ---°f and 
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with pi > 0, Yli=oPi = 1 an( i Vd > 0. We set w = Yli^Pi = /'(I) an d 7 := logci/logu;: we 
assume the super-criticality condition w > 1, so that 1 is an unstable fixed point of / and the 
sequence {/n(exp(/i))} n= o,i,... increases to infinity. Actually it is straightforward to see that 
the increase is super-exponentially fast and just a little more work leads to the existence of 
the free energy density of the model 

F(h) := lim -^log/„(exp(/i)). (1.2) 

We refer to [12} Appendix] for the existence of this limit as well as for various properties 
of F, such as the fact that F is a convex non-decreasing function. As a matter of fact, it is 
immediate to see that F(h) = for h < 0, while f(H) > for h > 0. The origin is therefore 
necessarily a critical point and it is actually associated to a localization transition (see [91 [12] 
and references therein for details on the statistical mechanics context, about which we are 
particularly concise here). 

One is then particularly interested in the free energy critical behavior, that is how f(/i) 
behaves near 0, which, in this case, of course reduces to considering h \, 0. What is argued 
in the aforementioned physical literature is that the expected critical behavior must be of the 
form 

F(h) WtfAQogh), (1.3) 

where the amplitude A : R — > (0, oo) is log w-periodic, which is of course compatible with 
^4(-) being constant. 

Remark 1.1. We point out that the hierarchical model in [9] has d = 2, po = (B — 1)/B, 
pi = and p2 = 1/B, with B G (1, 2) (to be precise in [9] the duaJ model with h replaced by 
log(-B — l) + h and B £ (2, oo) is considered, see [12] for the equivalence of these two models). 
The restriction to d = 2 is just for the sake of simplicity, but the choice of p\ =0 reflects 
a symmetry of the model that has no particular impact on the issue we tackle in this note: 
[9j [12] deal with disordered hierarchical pinning and in the presence of disorder the fact that 
pi > leads to a new phenomenon (but this is not not the case in absence of disorder). 

In [7] and with a substantially more detailed analysis in |10j the authors tackle the issue 
of establishing whether A(-) is trivial or not and of understanding the origin and size of the 
fluctuations. In |10| an argument is presented, in the Ising and Potts models framework, that 
is expected to capture the size of the critical amplitude oscillations, but it yields an amplitude 
that is smaller than the true one, obtained numerically, by several orders of magnitude. The 
authors argue heuristically about the reason of such a mismatch between their arguments and 
numerics and they point out the role of the geometry of the Julia set of the map associated 
to this model. In this direction they actually provide a relatively precise estimate of the 
oscillation amplitude by exploiting empirical estimates on the Julia set, even if they admit 
that the precision of their computation is surprising and presumably merely to be considered 
as a lucky circumstance \10\ p. 123]. Clarifying the relation between Julia set and oscillations 
is one of the main aims of this note. 

An important observation at this point is that there is a probabilistic representation of the 
partition function of the hierarchical pinning model [13] . Consider in fact a random branch- 
ing process, or Galton- Watson process, starting from one individual at time zero and with 
offspring distribution determined by the pj's weights. Therefore if one focuses on the number 
of individuals W n that are present at time n, then, conditionally on W n = k, W n+ \ is the 
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sum of k independent and identically distributed random variables for which the probability 
of being equal to j is pj . One directly verifies that 

f n (exp(h)) = E [exp(hW n )] . (1.4) 

1.2. Branching processes and the Harris random variable. It is a classical result that 
W n /w n converges almost surely to a non-degenerate limit random variable W (with a mass 
at zero if po > 0) . Moreover for every s£R 

lim E [exp(sW n /w n )] = E \exp(sW)] =: Ms) , (1.5) 

n— too 

and ip{-) extends to the whole complex plane as an entire function. In his seminal paper |16| . 
T. E. Harris pointed out (among several other facts) that 

logV>( s ) S ^°° s 7 L(logs) , (1.6) 

where L : R — > (0, oo) - the Harris function - is continuous and log u>-periodic. Harris was 
unable to show that L(-) is not constant, even if he was able to compute numerically the value 
of L(s) for d = 2, pi = 0.4 and pi = 0.6 up to six decimal digits, a remarkable achievement 
considering the date at which the paper was published. Later it became clear that L(-) 
does oscillate and that the amplitude of the oscillation is extremely small with respect to its 
average value (see in particular [3], but also [H O HZ]). A full understanding of this near- 
constant behavior is however still elusive. Actually analogous phenomena were recorded also 
in other mathematical fields like combinatorial enumeration, spectral properties of transition 
operators on fractals and more (see e.g. [201 EH ED E2] ) . Ultimately, this is not surprising 
because all these problems boil down to studying iterations of a map: in the Galton- Watson 
case for example (|1.4p tells us that the generating function of the law of W n is precisely f n . 

Certainly it has not escaped the reader that the qualitative properties of L(-) coincide with 
the (expected) qualitative properties of the critical amplitude A(-). However the quantitative 
connection between L(-) and A(-), beyond the common period w, is a priori not clear. In 
statistical mechanics terms, the Harris function L(-) emerges from the limit of the partition 
function in a particular vanishing limit of the pinning parameter s/w n : n — > oo and then 
s — > oo. A(-) comes also out of a limit of vanishing pinning parameter, but in this case 
n — t- oo is taken at fixed pinning parameter h > 0, only the leading Laplace asymptotic term 
is kept, and then h is sent to zero. 

Nevertheless, L(-) and A(-) coincide: 

Proposition 1.2. Given /(•) as above, the asymptotic relation (11.3P holds with A(-) a logw- 
periodic analytic function. Moreover L(-) = A(-) - from now on, they will be denoted by fl(-) 
- and, if we denote c n = J Q exp(inx)£l(x log w) dx the Fourier coefficients of£l{- log w), there 
exists a positive constant e such that \c n \ < exp(-en) for n sufficiently large. 

The notation we use does not highlight the dependence of fi(-) on /(•), but we stress that 
Proposition 11.21 says that given / one obtains a function f2. 

Proposition 11.21 can be proven in a rather direct way by exploiting a number of relations 
that one can find in the large literature devoted to the subject, but we have been unable to 
find the statement in this literature. The proof is in Section [2] but let us anticipate one of 
the main tools that is going to be of help for the sequel of the introduction: for every x > 

. oo 

F( l ogx) = ^iP£ + i ogx + Y^d'^Q(Mx)), (1.7) 



1 



OVIDIU COSTIN AND GIAMBATTISTA GIACOMIN 



with Q(x) := log(/(x)x d /pd)- This follows by observing that, if Fo(y) = y and F n (logx) := 
d- n logf n (x) forn = 1,2,... we have F n (log x) = d~ n logp d +d~ n Q(f n (x)) + F n - 1 (logx) [E]. 

To go beyond Proposition 11.21 we restrict to the d = 2 case. This restriction is made 
because we want to exploit directly the results in [UEJE], that develop only the case d = 2. 
It is certainly possible to generalize these works, but this would not add much to the purpose 
of this note at the expense of rather lengthy arguments. 

1.3. Julia set, Bottcher function and oscillations. If we set B{x) := exp(F(logx)) we 
directly verify that B solves the Bottcher equation B(f(x)) = B(x) d . This can be seen 
directly from (11. 7j) . that is B conjugates / and the monomial map x H > x d . B is increasing on 
[l,oo) and it is real analytic on (l,oo) (see the proof of Proposition II. 2\i . We can therefore 
set A := B~ l and A is analytic too. 

The central point of our approach is the following theorem, in which we use the notation 
D r := {z G C : \z\ < r}. In what follows log(-) is the natural logarithm with the choice of 
the negative semi-axis as branch cut: log(r exp(i#)) = logr + i0, for r > and 6 G (— vr,7r]. 
Let us recall [18] that the filled Julia set K = K(f) of the polynomial map / is the set 
{z G C : sup n |/ n (-2)| < oo} and that the Julia set J = J(f) is the boundary of the filled 
Julia set. 

It is immediate to see that J = {z G C : \z\ = 1} if f(x) = x d , that is if pd = 1 (recall 
(jl.ip ). but J is substantially more complicated if pd < 1. However note also that, for d = 2, 
/ is conjugated to the map g(x) = x 2 + c, c = (1 — (p2 — po) 2 )/4, via the affine transformation 
h(x) = x/p2 + (1 — Po — P2)/2p2, that is g = h^ 1 o / o h and that c G (0, 1/4): this implies 
that the Julia set is a simple closed curve (see for example [HI Ch. 3, Prop. 6.2]). 

Set d = 2 and < p 2 < 1. 

Theorem 1.3. The function A extends analytically to C \ Bi and to dO\ as a continuous 
function. Moreover there exist an entire function g, with g(0) = 1 and g'(0) > 0, and a 
\og2-periodic function u(-), which is non-trivial (i.e. non- constant), positive and analytic on 
the strip {z : |9(z)| < tt/2}, such that for \z\ > 1 we have 

A(z) = 5 ((logz) 1 /^(log(logz))) . (1.8) 

In addition: 

(1) the Julia set of f coincides with {A(exp(i7rt) : t G (—1, 1]}. 

(2) the function x i— > x 1 ^oj(log(x)), from the positive semi-axis to itself is invertible 
and its inverse is x ^ x 7 a(log(x)), where a is log w -periodic and fl(-) = c 7 a(- + 
log(c)), with c := P2/(2 — w). Moreover a(-), and therefore is analytic in the 
strip {z : |9(-z)| < ir/(2y)}, which implies that Jq° sw exp(27rmx/ log w)Q(x) dx = 
0(exp(— mT 2 /j')), for every 7' > 7 and n — > 00. 

To keep the statement simple the explicit series expansion for g(-) is postponed to Section [3l 
starting from fj3.4f) . In Section [3] one can also find an explicit construction for w(-). What 
we want to emphasize with Theorem 11.31 is the quantitative and explicit relation between 
oscillations and geometry of the Julia set. 

Remark 1.4. A more general but very implicit relation between Julia set and oscillations 
can be established by using the notion of Green's function of the monic polynomial map 
/ |18l P. 100-101], where monic means that f(x) = x d 4- 0(x d ~ 1 ). The Green's function 
u{z) := lim n _ >00 d~ n log |/ n (-^)| is well defined and in C°(C; M). It is actually harmonic ((d 2 + 
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dy )u(x + iy) = 0) except on the Julia set J(f) and it is of course identically zero on the filled 
Julia set K(f). Moreover u can be alternatively defined as the unique continuous function 
that vanishes on K(f), which is harmonic in J(/)^ and such that lim z _ 5 , 00 u(z) — log z = 0. 
Note the remarkable fact that u(-) depends on /(•) only through </(/). If now we observe 
that / is conjugated to a monic / via the linear map z \— > cz, c := p d *\ we readily see 

that J is cJ(f) and F(logx) = u(x/c) for x 6 (0, oo). Therefore the free energy is determined 
as unique the solution of a Dirichlet problem once J(f) (and the blowing up factor c) are 
given. As a consequence, the oscillatory function fi(-) is directly connected to the Julia set 
via the solution of this Dirichlet problem. 

Theorem 11.31 asserts that O(-) is not a constant: estimating the size of the oscillations is 
a challenging task. However, one does have explicit characterizations of uj and Vt that can 
be exploited to get precise numerical and also computer-assisted estimates on the Fourier 
coefficients of ui, with explicit error bounds. We address this issue in Section [3l 




Figure 1. The Julia set of f(x) = 1/4 + 3x 2 /4 plotted by using J = {A(exp(i-Kt) : t e 
( — 1, 1]} and (|1.8[) . after having determined A(-) as explained in Remark 13.11 a priori the 
method determines ^l(-) with precision in a neighborhood of 1 in the complex plane, as it 
is confirmed by the figure on the right, but the precision is quite satisfactory also far from 
1. We have superposed to this curve the points obtained by iterating (seven times) the 
pre-image map / , starting from the unstable fixed point 1. 

2. Proofs 

2.1. Proof is Proposition [TT2l We have the following representation of the Harris function 
[161 M L (-) given in ([LSI): for s > 

L(logs) = s~ 7 F(logV(s)). (2.1) 
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Figure 2. On the left the plot of fi(-) — fi, with SI := f ° gm Q(x) dx obtained by using 
f(x) = 1/4 + 3x 2 /4, that is the same map as in Figure [1] In particular the period logw is 
0.40546 . . . and we have Ti = 1.33381 . . . and the oscillation is 8.86 . . . x 10~ 8 . On the right 
we plot £!(■) — £1 — cisin(27r(- — xi)/logui), with ci the first Fourier coefficient and x\ the 
associated shift. This second plot makes clear the nearly sinusoidal character of £!(■) as well 
as the fact that the next Fourier coefficient is even much smaller: in fact the plot on the 
right is again nearly sinusoidal. 



It is worth recalling that (|2.ip follows by using the Poincare relation ip(ws) = f(^(s)), which 
is an immediate consequence of (|1.5p . and of F(log/(s)) = dF(log(s)) which is equivalent to 
the Bottcher functional relation (cf. § 1 1 . 3[) . In fact these two relations imply that if L(-) is 
the right-hand side of (|2.ip then we have 

L(logM) = F(lQgW " g))) = F(lQg/Ws))) = F(1 ° g(s)) = LQogs). (2.2) 

But (II. 7h directly implies that f(/i) ~ h and so 

~ F(log(^f(«; n «))) log(^(tO n a)) n->oo Tl . . . v rn s 

L{\ogs) = ~ L(n log w + logs = L(logs), (2.3) 

lOg(l/!(l0 n S)) (W n sp 

where in taking the limit we have also used (jl.6p . Therefore (12. ip is proven. 

Going back to the main argument, ifi : [0, oo) — > [1, oo) is increasing and ^(s) = l + s+o(s), 
for s small, so ip -1 (exp(/i)) ~ /i for /i small. Therefore, since L(-) is continuous (and periodic), 
we obtain from (|2.ip 

f(/i) ^%TL(log/i), (2.4) 

which is CL3]) with ^(-) = L(-). 

We are therefore left with the regularity properties of A(-) = L(-) which we now call f^(-). 
For this we first claim that that x i— > F(log(x)) extends to an analytic function in a cone 
C s := {z : |arg(z - 1)| < 6 and M(z - 1) > 0}, for a 5 > 0. 

The claim follows from by (|1.7j) because one can find c > 1 such that both | (f(z)z~ d /p^) — 
1| < 1/2 and 1/(^)1 > c for |z| > c. Therefore no singularity comes from the series in f)1.7|) if 
x G C,5, if (5 is sufficiently small, because by elementary estimates we see that the argument 
of f n ( x ) remains bounded by a constant smaller than ir for n < min{n : |/ ra (x)| > c}. 

Let us now look at the formula for fi(-) in (|2.ip with 1 + s G C5/ for < 5' < 5(< it): 
by periodicity it suffices to consider K(s) small, which guarantees that the (entire) function 
ip(s) = 1 + s + 0(s 2 ) takes values in the cone Cs and therefore Q(-) is analytic in the strip 
{z : \^s(z)\ < 5}, which directly implies that the Fourier coefficients of are smaller 
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than exp(27r<5'/ log w), for every 5' < 6 and n sufficiently large. This completes the proof of 
Proposition 11.21 □ 

The argument we have just presented effectively uses the analyticity of F(log(z)) for z in a 
truncated cone, that is for z G C$ and $t(z— 1) small. In order to get a better (in fact, optimal) 
estimate on the decay of the Fourier coefficients a precise knowledge of the Julia set is needed, 
in the sense that it is necessary for the truncated cone to be in the complement of the filled 
Julia set. This of course guarantees that f n (z) — > oo, but this a priori is not sufficient because 
one has to ensure also that the absolute value of the argument of f(z n )z~ d /pd, z n := f n (z), 
does not go beyond tt for every z = zq in the truncated cone (recall (II. 7D ). We therefore 
restrict to d = 2 and attack the problem from a somewhat different angle. 

2.2. Proof of Theorem 11.31 For the map / we consider oo is a super-attracting fixed point 
(see \18\ p. 44-45]). It is practical and customary to map the fixed point at infinity to a fixed 
point at zero by conjugation. Let us make this explicit for d = 2, so that f(x) = po+p\x+p2X 2 
and w = pi + 2p2 £ (1, 2). The fixed points of / are 1 (unstable) and P0/P2 < 1 (stable). The 
affine transformation x i— > l{x) 

l(x):=^- X ^, (2.5) 
2 — w 2 — w 

sends /(•) into the standard logistic map: 

/(/(r x (z))) = Xz(l-z), A := 2-w € (0,1). (2.6) 

The fixed points are mapped to 1 — 1/A (unstable) and (stable). We then map z to —1/z, 
so in the end we use the transformation q(x) := — l/l(x) (q : R \ {P0/P2} — > R \ {0}) to get 
to f (y) = (/(/(g _1 (?/))) and the new iteration 

y n+1 = f(y n ) = fl (2.7) 
A(l + y n ) 

and the unstable fixed point is now A/(l — A), while is stable (in fact, super-attracting). 
For later use, it is of help to note that q(l) = A/(l — A) and that 

8(h) := Y~\ ~ ^( ex P( /l )) ^ c a/i, with c A := ^^^2 ■ ( 2 ' 8 ) 

We are going to use the fact that there exists a unique function G, from Dj to the interior 
of K(f ), which is analytic and invertible, satisfying G(0) = and G'(0) = A and 

G _1 of oG(z) = z 2 . (2.9) 

Actually, the existence and uniqueness of such a map in a small disk around the origin is a 
general result (Bottcher Theorem [18, § 9]). G -1 is usually called Bottcher function for f and 
it is uniquely determined once we require it to be analytic near zero, with G _1 (0) = and 
(G _1 )'(0) = 1/A (this is a classical result: see [5] and references therein). The fact that G can 
be extended to the whole unit disk depends on the details of the map and, in our case, on 
the fact that A G (0, 1) [5J p. 1312]. 

The argument that follows is based on [5l Theorem 1] which gives an expression for </?(•) := 
1/G(-): there exists g(-) entire, with g(0) = and g'(0) = 1, and a non-trivial (non-constant!) 
periodic function w(-) of period log 2 (analytic on the strip {z : |9(z)| < vr/2}) such that 

<p{z) = ^ + g ((- log z^MM- log *))) , (2.10) 
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and we recall that 7 = log 2/ log w > 1. The properties of G(-) directly imply that ip(-) is 
analytic in the punctured unit disc and that it is a conformal (and invertible) map. Actually 
<p(-) extends by continuity to <9Bi [5, p. 1313]. Moreover g(-) is invertible in the range of its 
argument in (I2.10p . 

Let us point out that (I2.10p implies (|1.8p . For this let us us go back to the (|1.7p . with 
which we have defined B and then A, see § 11.31 As we have already remarked in the proof of 
Proposition 11.21 l/n(-z)| > \z\ for \z\ sufficiently large and this directly implies the analyticity 
of B{) in a neighborhood of infinity, as well as the existence of it inverse A(-) in the same 
neighborhood, since B(z) ~ P2Z for \z\ —> 00 (we remark that B(-) coincides with the function 
ip in |181 p. 98]). Going (backward) through the conjugation that we have performed we see 
that \/B(q~ 1 (z)) ~ z/X for \z\ large and that z \-t 1/B(q~ 1 (z)) solves the Bottcher equation 
(in the sense that it has the same property as G _1 in (|2.9|) ) . It must therefore coincide 
with G -1 and this extends the domain of analyticity, in fact bi-analyticity, of A(-). Again, 
by going carefully through the backward conjugation we see that g(-) = 1 + (X/p2)g(-), so 
(|l,8p is proven. The fact that A(-) extends as a continuous function to <9Bi follows from the 
analogous property for <£>(•). 

Point (1) of Theorem 11.31 follows from the general theory |18t Theorem 9.5]. In fact the 
Bottcher function £>(•) maps (in a bi- holomorphic fashion) K(ff to B? and A{dB{) = J(f). 

Let us move to point (2) and let us start by remarking that 

1 1 

F(h) = lim — log/„(exp(/i)) = - lim — log f „ (q{exp(h))) . (2.11) 

Therefore, in view of (|2.8p we need to control f n ((A/(l — A)) — 5) for n — > 00 and then 5 
small. But (j2.9|) tells us that for every n 

f n (G(x)) = G(x 2n ) . (2.12) 

and therefore 

'•(iV'HHiV'f)- (2i3) 

It is practical to set (in particular for S > 0) 

i{5) := -logG- 1 ^^-^ . (2.14) 

We claim the following: 



1-A 



Lemma 2.1. The map x 1— > x l I ^uj (log x) is a bijection from (0, 00) to itself and its inverse 
can be written as x \-> x 7 a(log(a;)), with a{-) log w -periodic. Moreover a(-) is analytic in the 
strip {z : Q(z) < ir/(2j)}. Finally for c = ((1 — A)/A) 2 we have 

l{5) S ^>° {5cpa (log(<5c)) . (2.15) 

Proof. We use the properties of cp(-) (or G(-)) and (|2.10p . Therefore we see that if we set 
z = exp(—£(5)) we have 

<p (exp {-1(8))) = x _\~^ x)d =■■ ^ + v(S) , (2.16) 

and v(6) ~ c6, c given in the statement. Thus we have to solve 

£(5)^u; (log(£(6))) = g- 1 (v(5)) ( ^° cS) , (2.17) 
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at least for 6 small. We see therefore that it is a matter of inverting x h-> x 1//7 cj(log x) =: /3(x). 
First, observe that /3 : (0, oo) — >■ (0, oo) is invertible, because ip is invertible. Moreover /?(•) 
is (real) analytic, hence /3 _1 (-) is too. Then set a(logy) := /3 _1 (y)/y 7 , so that a(-) is real 
analytic and 

a [ -\ogx + \ogu{\ogx) J = — r r. (2.18) 

V7 / wT(logx) 

From (I2.18P we directly see that, since u(-) is log 2-periodic, a(-) is (7 -1 log 2)-periodic (that, 
is log w-periodic) and, since u(-) is not constant, a(-) is not a constant either . Therefore the 
statement is proven, except for the domain of analyticity. To this end note that the argument 
above was given by restricting to the real axis, but we do know that uj(-) is analytic on the 
symmetric strip of half-width tt/2 and therefore /?(•) is analytic in the positive half-plane 
{z : 3i(z) > 0}: all we need to know is the fact that it is invertible in this domain. But 
since <p(-) is defined and invertible in the punctured unit disk, z 1— > y?(exp(— z)) is defined and 
invertible in {z : di(z) > 0} and, thanks to ()2. 16j) and (|2.17p . we see that also /?(•) is invertible, 
at least if we restrict to the truncated cone of the points z such that |Arg(z)| < vr/(l + e), 
for any choice of e > 0, and 5ft(z) > sufficiently small. Once again, since a(-) is periodic, 
analyticity on {z : \^s(z)\ < 71-/(27(1 + e)) and ^St(z) < —C}, for s > arbitrary and C = C(e) 
large, implies analyticity on the whole strip {z : \^s(z)\ < 77/(27)}. □ 

Let us complete the proof of (2) and for this let us go back to ()2.11|) and use (|2.8p and 
([XT3D to see that 

F <' ! > = -J™ F losG ( G " (t^a-' 5 *'") 2 ) (2 - 19) 

A direct consequence of Lemma [2. II and of Hm^x $(h) = is that G _1 ((A/1 — A) — 6(h)) < 1 
for h sufficiently small and therefore, since G(0) = and G'(0) > 0, for such values of h, we 
have 

F(h) = - Hm 1 log G- 1 - 6(h)) 2 = - log G- 1 - 6(h)) = l(S(h)) ■ 

(2.20) 

We now recall once again (|2.8p and apply Lemma 12.11 to obtain 

F(h) h >t (c x chya(log(c x ch)) , (2.21) 
and we compute c\c = P2/A = c, from which we identify O(-) in terms of a(-) and c. □ 



Remark 2.2. Of course one can upgrade the proof we just completed to include Proposi- 
tion [L2l that is to deal with the asymptotic behavior of V'(')- This is straightforward, albeit 
a bit lengthy: we sketch the main steps and make some comments. First of all we have 

iP(s) = lim q- 1 (f n (q(exp(s/w n )))) , (2.22) 

n— >oo 

and we notice the analogy with (|2. 13|) and we write 

f n (q(exp(s/w n ))) = G U' 1 L^- - 6 n \ j , (2.23) 



where 



5(e^ -1) rt ~°°CA— . (2.24) 
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It is now a matter of applying Lemma 12.11 the net result is the following representation for 
the generating function of the Harris random variable: 

= q- 1 (G (exp (-(scf a(log( S c))))) , (2.25) 

where we recall that c = P2/A. Since g _1 o G(x) for x \ behaves like a constant times l/x, 
one readily recovers (jl.6p (let us remark that in fact (|2.25p is equivalent to (|2.ip ), 

3. On quantitative estimates 

3.1. Estimating A detailed numerical approach to f2(-) can be found in [3] for f(x) = 
Ax + (1-A)x 2 , for A = 0.1, 0.2, . . . , 0.9, but without explicit error bounds. The (two) methods 
they employ however do allow an explicit control of the error, with a non-trivial amount of 
work. Notably, by using what they call Bottcher method, which is based on (|2.ip . one uses 
(jl.7p for f(-), for which it is straightforward to control the error when one truncates the series, 
and one can easily set up an iterative procedure to get the Taylor coefficients of ?/>(•) with a 
control on the remainder by fixed point arguments. Note that it is sufficient to control the 
error for s £ [c, cw], for a conveniently chosen c. Similar ideas are developed in detail just 
below for w(-), on which we focus, and we recover Q(-) from oj(-). 

3.2. Estimating u>(-). By recalling (I2.10p and by the properties stated right after that for- 
mula we see that 

a;(log S ) = s-V 1 (V(exp(-s)) - ^) , (3.1) 

where (see [5]) 

• g(-) is the unique solution to 

g(y) = A(g(y/w)), with A(y) := wy + Ay 2 , g(0) = 0, g'(0) = 1; (3.2) 

• ip{-) is the unique solution of 

ip(y 2 ) = P(tp(y)) , with P(x) = \x{l + x) , (3.3) 

under the condition that (/?(•) is analytic in the punctured unit disc and that (p(z) ~ 
l/(\z) for z going to zero. 

We recall the relation 

<?(•) = l + yg(-), (3-4) 

and <?(•) is the function appearing in Theorem 11.31 

In principle, it is now a matter of exploiting (|3.2|) and (|3.3|) to get enough terms in the power 
series of g and ip and of controlling the remainder, to obtain a sufficiently good approximation 
of oj on [x, x + log 2], for some x, from which we can extract the Fourier coefficients of u by 
controlled numerical integration. While in principle this whole procedure is straightforward, 
in practice it is quite non-trivial given the fast decay of the coefficients and the fact that even 
the first coefficients are extremely small. 

Remark 3.1. Once uj(-) is precisely estimated over an interval of length log 2 - the period - 
it is of course known with the same precision over 1R. The plot of the Julia set in[T]is obtained 
by following the same principle: we will not perform explicit estimates for uj(-) on C and we 
content ourselves with remarking the somewhat surprising precision of such a procedure, see 
Fig. [T] and its caption. The graph in Fig. [T] has been obtained by keeping 250 terms in both 
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the series for g and ip, and by performing three times the backward iteration procedure that 
we explain just below. 

3.3. Approximating g and g _1 . From (|3.2p we see that the power series for g has positive 
coefficients. One can obtain the first n coefficients by setting g^(x) := x + g^x 1 + . . . + g n x n , 
by computing the polynomial g^(wx) — A(gr n )(x)) and by setting to zero the coefficients of 
the terms of degree smaller than n + 1. This determines g%, . . . ,g n . Then we define q via 
g(x) = g( n ){x) + y n+l q(x) and use (|3.2p to obtain 

q(x) = Q{x) +p(x)q(x/w) + 2( ^ +1) x n+1 {q(x/w)) 2 , (3.5) 

where Q is a polynomial of degree n — 1 and Q(0) ^ 0, p(x) := w~ n + 2\w~ n - 1 g {n) (x/w). 
With the notation \\q\\ £ '■= max^^e] \q(x)\ we directly obtain that if we can exhibit a > 1 
and e > such that 

c e := 1 - ||p|| e - a -^£ n+1 \Q(0)\ > and ||Q|| e < ac £ \Q(0)\ , (3.6) 

then \q(x)\ < a\Q(0)\ for |x| < e. Therefore, by the positivity of the series coefficients, we 
have 

g {n) {x) < g(x) < 9{n) (x) + a|Q(0)|x n+1 , (3.7) 

for x S [0, e]. By inverting the two polynomials that bound g from below and above, and 
by taking the Taylor expansion to order n of these two expressions one directly recovers the 
power series for g _1 truncated at n, with an explicit control on the rest. 

A performing way to improve this approximation of g" 1 is the following: from (|3.2p we 
obtain 

&-\y) = wg-^A-^y)) , (3.8) 

where 

A-Hv)= (3 . 9) 

A^ 1 is concave, with slope 1/w at the origin and from this we get that if g -1 (y) = p{y)+fo{y), 
with p(-) > and ro a remainder like above, we get 

\g-\y)-wp[A- 1 (y))\ < ^IroCA- 1 ^))! < Cw (^)* . (3.10) 

Note that the new remainder is better both because it improves by a factor w 1 ~ k the estimate 
in the interval in which we have the estimate for ro (with no a priori condition on the argument 
of ro) and because it yields an explicit estimate on an interval that is w times larger. We 
sum up this argument: 

Fact 3.2. If g^ 1 =p + tq, with p{x) > 0, |ro(x)| < Cx k for x e [0, xq] and if we set 

r n (y) ■= g-\y) - w n p ((A~ 1 ) on (y)) , (3.11) 

then \r n (y)\ < w n \r ((A' 1 ) on (y))\ < Cy k /w^ k ~ 1 ^ n , for y <E [0, w n xo] (actually, even on a 
much larger set, since A on (y) becomes much larger than yw n for n large). 
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3.4. Approximating <p. One starts by guessing the (Laurent) series coefficients for tp, by 
using (13.3|) . and it is not difficult to see that 

= ^ " \ + P(x) + y 2n+1 r(y) , (3.12) 

where r is analytic and p{x) is a polynomial of degree 2n — 1 that contains only odd powers 
of x. Again by using (13 .3|) we extract an equation for r: 

r(x)pi(x) + p (x) - \x 2n+2 r n {xf + r n (x 2 )x 2n+2 = , (3.13) 

where p\(x) is a polynomial of degree 2n, containing only even powers, with pi(0) = —2; 
Po(x) contains also only even powers and it is a polynomial of degree 2n — 2. Then one can 
show that if one can exhibit a > 1 and e > such that 



C e := 1 



1 


e 2n+2 (a\ 


Po(0) 


Pi 




Pi(0) 



+ 1 > and 



Po 
Pi 



< aC F 



Po(0) 



Pi(0) 



then 



</?(x) 



1 



+p(x; 



< ^bo(0)|/" +1 



(3.14) 
(3.15) 



1 

M 2 

for (0,e]. 

The approximation of <p can be greatly improved by backward iteration, as for g _1 : from 
(1331) it follows that 



<p(y) = Qi^iy 2 )) , where Q(x) 



i + ¥ -i 



(3.16) 



Note that Q(-) is concave, Q'(0) = 1/A and Q(x) = x implies x = or x = (1 — A)/A. So if 
ip = (p + ro, with tp(x) > and |ro(x)| < Cx k for x G (0,xq], then 



1 



C 



2A- 



(3.17) 



for y E (0, y^xo]. This can also be seen by using max^o |Q'(x)| = 1/A. Note that (|3.17[) 
implies that Q(^p{y 2 )) is better that <p(y) in approximating tp(y) in the sense that the new 
remainder (r±(y) := ip(y) — Q(<p(y 2 ))) is smaller than the previous one if y is in the original 
interval, that is y E (0, xo], if x o < A 1//fc , and one has an explicit estimate for r\ in a larger 
interval. To sum up: 

Fact 3.3. // ip = <p> + tq, with tp(x) > and |ro(x)| < Cx k for x S (0, xo], then if we set 

r n (y):=^y)-Q on &{y 2n )) , (3.18) 
we have \r n (y)\ < ^|r (y 2n )| < Cy 2nk /\ n , for y G [0,xj /2 "]. 



3.5. Example of computation. By applying the arguments of the previous subsection for 
A = 2 - w = 1/5 we have 

, s (2-w)x 2 2(2-w) 2 x 3 (2-wf(h + w)x i 

g(x) = X + ) ±- + \.. ' ; . + 



(to — (u> — l) 2 it; 2 (l + w) (w — l) 3 w 3 (l + w)(l + w + w 2 

2(w - 2) 4 (7 + w(3 + 2w))x 5 



{w - l) 4 w 4 (l + w) 2 (l + w 2 )(l + w + w 2 ) 



+ r 5 (x), (3.19) 
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with |rs(x)| < 5|x| 6 /10 6 for \x\ < 2.5. From this we extract 



g 1 (x) = x 



Xx 2 



+ 



2X 2 x 



2^3 



A 3 (l + 5w 2 ) x A 



{—l + w)w (— 1 + w) 2 w(l + w) (— 1 + w) 3 w 2 (l + w) (1 + w + w 2 ) 

2A 4 (3 + 2w + 7u> 3 ) 



+ 



m3 I 



(-1 + w) A w{\ + w) 2 (1 + w 2 ) (1 + w + w 2 ) 



+ r 5 -(x), (3.20) 



with |r 5 {x)\ < 5\x\ °/10 for |x| < 2. Moreover 
<Pi?) = ^-^ + ^(-2 + A)x + T ^(-16 + 4A + 4A 2 -A 3 )x 3 



+ 



where \R*(x)\ < x 7 /W for x £ (0,9/10]. 



-2 + A)A(l6-4A-4A 2 + A 3 ) 
1024 



+ R 6 (x), (3.21) 




Figure 3. The plots of approximations of w(-) — ZJ, cf. p.22p . for A = 1/5 are obtained 
by using the expansion of g _1 to 5 th order and <p to 6 th order, and by applying five times, 
on the left, and seven times, on the right, the backward iteration procedure. The case on 
the right is used for the quantitative estimates, and the interval [—2 — log 2, —2], in which 
we have good error estimates, is chosen. 



Aiming at obtaining a good approximation of u; we choose s G [l/(2e 2 ), 1/e 2 ], so exp(— s) G 
[0.873 . . . , 0.934 . . .]. For this interval of values we obtain <p(exp(-s))-(l-A)/A G [0.48 . . . , 0.92 . . .] 
with an error of at most (0.93) 7x2ni /(10(0.2) ni , where rii is the number of backward iterations 
employed, cf. (I3.18p . Similarly, the error on g _1 in the range of values under consideration 
is uniformly bounded by (5/10 4 )(0.93) 6 /(1.8) 5 ™ 2 . Choosing n\ = 7 makes the error on ip 
smaller than 10 -23 and setting also ri2 = 7 leads to an error on uj that is uniformly bounded 
by 5 x 10 -13 . This allows to estimate the average value of to and the first Fourier coefficient 
with a precision of at least 10 -11 : uj(x) =to + gi sin(2-7r(x — Xq) / log 2) + . . . 

i r 2 

U:=- / oo{y)dy rj 4.45140273002 (3.22) 

log 2 y_ 2 -log(2) 

and g x w 5.938 x 10~ 8 . 

At this point if one wants to recover fi(-) one has to perform the inversion step in Theo- 
rem [L3j2). It is not difficult to realize by considering (12.180 that, given the small size of the 
oscillations, a good approximation of a{y) (defined in Theorem ll.3( 2)) is Xjuf'^y — 7 log a;)). 
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Of course a control of the error requires an attentive (but elementary) analysis. By perform- 
ing explicitly the case under consideration and choosing for example p\ = 0, P2 = w/2 and 
p = 1 - P2 we obtain H = 1.01288677326 and the first Fourier coefficient is 1.59 x 1(T 8 . 
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